February 6, 2008 23:8 



WSPC - Proceedings Trim Size: 9.75in x 6.5in qscpll'suraud 



Size and charge effects on the deposition of Na on Ar 

P. M. Dinh 1 , F. Fehrer 2 , P.-G. Reinhard 2 , and E. Suraud 1 

1 Laboratoire de Physique Theorique, UMR 5152, Universite Paul Sabatier, 118 route de 
Narbonne, F-31062 Toulouse Cedex, France 

2 Institut fur Theoretische Physik, Universitdt Erlangen Staudtstr. 7, D-91058 Erlangen, 

Germany 

We discuss the dynamical deposition of the Na atom, the Na+ in and the Na— 6 cluster 
on finite Ar clusters mocking up an infinite Ar surface. We analyze this scenario as a 
function of projectile initial kinetic energy and of the size the target cluster. 

1. Introduction 

Clusters on surfaces have motivated many studies over the past decades and still 
remain a topic of great interest, 1 especially in relation to the synthesis of nanostruc- 
tured surfaces. 2-5 Indeed, it turns out that it is possible to make a direct deposition 
of size selected clusters on a substrate. 6 ' 7 The deposition process may lead to a 
significant modification of the cluster, in terms of its electronic structure and ionic 
geometry. This is a consequence of the impact of the interface energy, the electronic 
band structure of the substrate, and the surface corrugation. These questions have 
already been widely investigated in great detail, especially from the structural point 
of view and both from the experimental 8-10 and theoretical 1 1-18 sides. The situa- 
tion is somewhat different concerning the deposition process itself. Its theoretical 
description remains mostly limited to molecular dynamics (MD) approaches, which 
implies that a proper description of electronic degrees of freedom is missing. The 
reason for this defect is simple. The presence of a substrate makes the experimen- 
tal handling of clusters easier but strongly complicates the theoretical description 
because of the huge number of degrees of freedom of the surface. The MD then 
provides the cheapest way to access, at least in a gross way, the dynamics of the 
substrate. It is nevertheless crucial to try to account for the surface's electronic 
degrees of freedom, for example when non adiabatic processes are involved. 

As a first step in the direction of a fully microscopic dynamical approach, one 
may consider relatively simple cluster/substrate combinations. This is for example 
the case of the deposit of a metal cluster on an insulator surface. The surface can 
then be included at a lower level of description, which simplifies the handling, as 
was, e.g., explored for the case of Na clusters on NaCl in. 18 ' 19 Cluster electrons 
were there described by means of Density Functional Theory and the coupling to 
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the surface was achieved via an effective interface potential, itself tuned to ab-initio 
calculations. 14 Such an approach implies a total freezing of the surface itself, which 
sets severe limitations on its applicability. A somewhat better description of surface 
degrees of freedom in a still limited/simplified way can be achieved by considering 
again a rather inert substrate, but allowing for a minimum of dynamical response 
of the substrate. Such a model was recently proposed for describing sodium clusters 
embedded in rare gases. 20-22 The method used was "hierarchical" , which is justified 
by the moderate interactions between cluster and surface. The interface can then 
be treated at a simpler level than the cluster's degrees of freedom (microscopic 
treatment of cluster and classical treatment of environment with explicit account of 
dynamical polarizability effects), in the spirit of the coupled quantum- mechanical 
with molecular-mechanical method (QM/MM) often used in bio-chemistry. 23-25 

Such hierarchical methods, although much simpler than a fully microscopic ap- 
proach, require sophisticated modeling and are thus restricted to finite systems. 
Nonetheless, the calculations in 21 ' 22 were carried forth to a sufficiently large range 
of sizes to see the appearance of generic behaviors on the way towards the bulk. In 
the complementing case of deposit on a surface, we also consider finite substrates, 
as model cases for a surface. This is an acceptable compromise once the impact of 
the finiteness of the substrate has been properly analyzed. For a first exploration, 
we went one step further and restricted the present analysis to the even simpler 
case of the deposition of a single Na atom on a finite Ar cluster, adding one test 
case in which we shall explore the case of a finite Na cluster. 

The goal of this paper is thus mostly the study of the dynamics of deposition of a 
sodium atom (projectile) on Ar clusters of various sizes Arjv(target, A = 43, 86). We 
shall analyze the behaviors of both the atom and the cluster, especially as a function 
of deposit velocity and also consider size and charge effects (deposit of a Nag cluster 
and a Na + ion instead of a Na atom). The paper is organized as follows. Section 
[2] gives a short presentation of the model used. A few more details on the model 
can be found in the appendix, where basic formulae and parameters are recalled. 
The following sections successively address the dependence on substrate size and 
on projectile velocity. We finally discuss the example of a true cluster deposit and 
of charge effects. 

2. Model 

The model was presented in detail in 26 and we provide the basic formulae in the 
appendix. We recall briefly the basic ingredients. The Na cluster is described in 
terms of time-dependent local-density approximation for the electrons coupled to 
molecular dynamics for the ions (TDLDA-MD) , a scheme which has been exten- 
sively validated for linear and non-linear dynamics in free metal clusters. 27,28 The 
electron-ion interaction is treated by means of a soft, local pseudo-potentials. 29 
Each Ar atom is described by two classical degrees-of-freedom: its center-of-mass 
and its electrical dipole moment. The explicit account of Ar dipoles allows to treat 
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the polarizability of the atoms dynamically, with help of polarization potentials. 30 
Atom-atom interactions are described by a standard Lennard- Jones potential. For 
the Ar-Na + interaction we employ effective potentials from the literature. 31 The 
electron- Ar core repulsion is modeled in the form proposed by, 32 with a slight final 
readjustment of the parameters to the NaAr molecule as benchmark (bond length, 
binding energy and optical excitation). A Van der Waals interaction is also added 
and computed via the variance of dipole operators. 20,26 ' 32 




Fig. 1. Atomic configuration of Arg6 with respect to the single Na (white ball). Six Ar atoms are 
emphasized by dark gray balls for the discussion on Na+ deposition in Arga (see text). 

As explained in the appendix, the starting quantity is the total energy, con- 
structed from the various pieces discussed above. The corresponding equations 
of motion can then be derived by standard variation, which leads to the time- 
dependent Kohn-Sham equations for the cluster electrons. One furthermore obtains 
Hamiltonian equations of motion for the classical degrees of freedom (Na + ions as 
well as Ar atom positions and dipoles). The initial condition is provided by the 
corresponding stationary solutions. 

The definition of the Ar-Na configuration requires some specific handling. Indeed 
the Na atom can be initially placed above either an Ar atom or an interstitial site 
of the surface layer. In the case of Ar43 (see Figure [2]), the first option has been 
chosen, and for Ar§6, the second initial configuration has been used, see Figure [TJ 
The Ar§6 is obtained from the stable free Ar§7 cluster where an outer Ar atom has 
been removed. The resulting Ars6 is then rotated in order to present a flat surface 
to the impinging Na atom, this way simulating the flat interface provided by an 
infinite surface. The starting configuration for the deposit of Nag on Ar43 is similar 
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to that of the single Na atom. This means the top Na ion (the Na6 is composed of 
a ring of 5 ions and an outer ion) is placed above an Ar atom. 

The numerical solution proceeds with standard methods as detailed in. 28 The 
Kohn-Sham equations for the cluster electrons are solved using real space grid tech- 
niques. The time propagation proceeds using a time-splitting method. The station- 
ary solution is attained by accelerated gradient iterations. We furthermore employ 
the cylindrically-averaged pseudo-potential scheme (CAPS) as introduced in, 33,34 
an approximation justified for the chosen test cases. In the following, the symmetry 
axis is denoted by the z axis. It should nevertheless be noted that the dynamics of 
the Na + ions as well as that of the Ar atoms are treated in full 3D. 



3. Dynamical deposition of Na on finite Ar clusters 

3.1. An example 

We first consider, as an example, the case of the deposition of a Na atom on an 
A1-43 cluster. This is illustrated in Figure [2] where a few snapshots of the deposition 
process are presented. In that test case, the Ar cluster presents to the Na atom a 
rather small surface area of only 5 atoms and the the Na atom is initially positioned 
above an Ar atom (which shows strong core repulsion). These two features differ 
as compared to the case of Ar 86 illustrated in Figure Q] and discussed later on in 
Section 13.31 They should a priori make an attachment of the Na sodium more 
difficult here. Still, the calculation shows that a faint binding takes place in this 
case within typically 3 ps. The Na atom loosely attaches to the Ar cluster surface, 
since it keeps on bouncing over 7 ps with a decreasing amplitude and around an 
average position of about 8 ao, even a bit closer than the NaAr dimer bond length 
of 9.5 ao- Mind that the initial Na kinetic energy Eq was 13.6 meV, in between the 
5 meV binding energy of the NaAr dimer and the 50 meV of the bonding in Ar 
bulk. Thus there is no surprise to observe the creation of a transient NaAr bond 
and a very slight rearrangement of the Ar cluster, which takes place over a time 
scale of order of 5 ps. Actually, the Na atom is accelerated during its fall and gets, 
before the hit, a kinetic energy almost three times higher than Eq (this depends on 
its initial separation with the Ar first layer). Then after the hit, the extra kinetic 
energy is absorbed by the Ar cluster both at the side of kinetic energy (about 15 
meV) and in terms of potential energy in its structure rearrangement. The residual 
kinetic energy of Ar atoms as well as the potential energy consumed in the Ar 
cluster rearrangement strongly depend on the available energy and the number of 
degrees of freedom. It is thus interesting to study the effect of variations of both 
these parameters. 

3.2. Dependence on kinetic energy 

We first consider the influence of the initial kinetic energy of the Na projectile. 
The results are displayed in Figure [3] for the deposition of a Na atom on Ar§6 at 
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Fig. 2. Example of dynamical deposition of a Na atom on an Ar43 cluster, with initial kinetic 
energy Eo = 13.6 meV. 

various energies (E = 13.6 to 870 meV). It should first be noted that, whatever 
the initial kinetic energy no electronic emission is observed. We shall thus restrict 
the discussion to ionic and atomic degrees of freedom. In the case of the two lowest 
impact energies, one observes some bounces, and finally a binding of the Na to the 
surface at a typical distance of about 8 — 10 ao from the first layer. However for 
higher Eo, the Na is reflected by the Ar cluster, and no attachment of the Na is 
observed, although the Na is initially located above an interstitial position. This 
result can be confirmed by computing the Born-Oppenheimer surface of Na on an 
infinite Ar surface. The ground state surface exhibits a faint minimum around 7 ao 
above the Ar surface, which is fully compatible with our dynamical analysis on the 
finite Ar 86 cluster. Closer to the Ar surface, the Born-Oppenheimer surface exhibits 
a strong repulsion reflecting the core repulsion between Na and Ar. When the impact 
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Fig. 3. Dynamical deposition of Na (thick lines) on Arg6 for four different initial kinetic energies. 
The z-coordinates (left) and the kinetic energies (right) are presented as a function of time. 

kinetic energy becomes too large the Na atom thus "misses" the faint minimum and 
directly explores the strongly repulsive part of the potential. The Na atom is then 
reflected by the cluster as observed in Figure The energy threshold for neutral 
Na sticking seems to lie between 0.05 and 0.2 eV. Note that in this range of kinetic 
energies, the Ar cluster is not affected very much. This is particularly visible on 
the time dependence of its kinetic energy, which exhibits a rather smooth energy 
transfer from the Na to the Ar cluster. We also observe the propagation of a soft 
shock wave in the substrate. This participates to the slight heating of the Ar cluster 
(to a few tens of K), except for the highest energy where some deep Ar atoms are 
emitted because of the stronger wave propagating through the layers. 

3.3. Dependence on Ar cluster size 

We now consider the influence of the number of degrees of freedom on the capacity 
to dissipate the available energy. In practice, this amounts to test the influence of the 
Ar cluster size for constant initial kinetic energy of the Na atom. The comparison 
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Fig. 4. Deposit of a single Na atom (Eo = 13.6) on Arjv systems of various sizes (TV = 43 and 
86). Are presented as a function of time, the z coordinates (left) and the kinetic energies (right) 
of the Na atom (thick lines) and the Ar cluster (thin curves). 



is presented in Figure [4] where we plot Na and Ar positions and kinetic energies 
as a function of time during the deposition process of a Na atom (with initial 
kinetic energy Eo of 13.6 meV) on Ar43 and Arg6- As already mentioned, two main 
differences between both cases are to be noted. First the Na atom is initially above 
an Ar atom of the Ar43, at a distance of 20 ao- In the case of Arg6, the Na atom 
is above an interstitial site and starts at a smaller distance from the Ar first layer, 
namely 15 ao. Due to their different sizes, the effective Ar surface exhibits only 5 
atoms in the case of Ar43 but 12 Ar atoms for Ar§6- 

Similar oscillating patterns are observed in both cases. 2 ps after the impact, 
almost all kinetic energy of the Na atom is transfered to the Ar cluster, while the 
Ar kinetic energy seems to reach an equilibrium and oscillates around a mean value 
slightly higher than the initial kinetic energy Eq . Mind that the transferred energy 
is distributed equally over all Ar atoms. This means that, as expected, the larger 
the target cluster size (thus the more available degrees of freedom) , the smaller the 
relative energy shared per Ar atom, and so the more moderate the perturbation at 
their side. Note, furthermore, that more Na kinetic energy is available at the time 
of impact (maximum of Na kinetic energy in the right panels) for the case of Ar§6- 
This means that the Na atom is accelerated faster in that case as compared with 
Ar43. The reason is that the larger cluster and its larger surface area provide more 
attraction from polarization potentials. Nonetheless, in both cases, the Na atom 
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loosely binds to the surface, at about the same distance of about 8 oo- 
4. Example of a finite cluster deposit 
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Fig. 5. Example of deposit of a finite Na6 cluster on Ar43. 

In order to complement the analysis in terms of Ar cluster size, it is also inter- 
esting to consider the case of the deposition of a full Na cluster, instead of a single 
atom. This again enhances the number of degrees of freedom. We shall consider 
comparable impact kinetic energy per Na atom. In Figure [5] we consider an exam- 
ple of a deposition of a finite Na cluster (Nae) on Ar43. The initial kinetic energy 
of the cluster is -Ekino = 13-6 meV per Na atom. As in the case of a single Na atom 
(Figure [2]) the pinning process proceeds stepwise with a slight bounce before the 
metal cluster finally attaches to the rare gas cluster. At variance with the case of 
a single atom, though, one can furthermore analyze the evolution of the shape of 
the deposited Na cluster. The considered Na6 cluster is primarily strongly oblate, 
consisting of a pentagon of 5 Na atoms topped by one central atom. One can see 
that during the deposition process, the cluster shape is little affected. This might 
have been expected in view of the ideal "flat" shape of the Nag cluster which already 
presents a large contact area to the Ar surface. But we also found for other, and less 
favorable, geometries (e.g. the nearly spherical Nas cluster) that the cluster shape 
remains basically intact during deposition. The details of this scenario of course 
depend on the size of the Ar cluster target and on the initial kinetic energy of the 
impinging Na cluster, but qualitatively the example displayed in Figure [5] turns out 
to be quite typical. These results will be presented elsewhere. 38 
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Fig. 6. Same as in Figure [3] in the case of Na + on Ars6- 



As a final point, we want now to analyze the influence of charge on the deposition 
process. We come back to the simple case of a single Na atom and consider the 
deposition of the corresponding Na + ion on Ar§6, at various initial kinetic energies. 
Born Oppenheimer calculations of Na + in contact to an Ar surface show that the 
attachment is much stronger (and closer to the surface) than in the case of the 
neutral species. This can be easily understood by remembering the key role played by 
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the attractive Ar polarization potentials. The finite charge of the Na + ion strongly 
polarizes the surface and thus enhances the binding as compared to the neutral 
case. Of course short range repulsion remains present but will take the lead only on 
shorter distances. The attachment is thus expected to be stronger and closer to the 
surface. 

Figure [6] confirms the expectation. The Na + ion is practically swallowed by the 
Ar cluster and the Ar cluster itself undergoes stronger rearrangements. For the 
smallest Eq presented in the top panels of Figure two of the six surface atoms are 
finally ejected from the Ar cluster (they are visible by the light lines going straight 
through Ar layers in the z-coordinate panel and having reached the lower end at 
7 ps) . The four other Ar atoms of the (strongly perturbed) surface remain bound 
to the whole cluster and participate to its strong rearrangement. In the case of 
intermediate Eq, the six Ar atoms are lifted from the first layer but still stick to 
the edges of the Ar cluster. The fact that the Na + goes deeper for Eq = 54.4 meV 
seems to be due to the final ejection of an Ar atom. In the latter three cases, these 
six atoms get about 10% of the total Ar kinetic energy. Finally, for the highest Eq, 
the first layer just explodes after the impact. The six outmost Ar atoms absorb up 
to 50% of the vertical shock in term of kinetic energy transfer, and then follow a 
radial motion in an horizontal plane (see the thin horizontal lines around z = in 
the bottom left panel of Fig|6|). In all cases, the removal of the six atoms from the 
first Ar layer, either outside or at the edges of the Ar cluster, allows some Ar atoms 
deeper inside the substrate to move upwards, thus leaving a large vacancy, so that 
the Na + can penetrate even under the second layer. These results suggest that the 
number of layers under the Ar surface probably does not play an important role in 
the Na + inclusion. More important is the surface's mobility for rearrangements. 

6. Conclusion 

In this paper, we have presented results on the deposition of a Na atom, a Na + ion, 
and a Nag cluster on a dynamically polarizable Ar substrate represented by finite 
Ar clusters of various sizes. We have used time-dependent density-functional theory 
for the Na electrons coupled to molecular dynamics for the treatment of Na ions 
and Ar atoms. We have presented systematic results as a function of Ar cluster size 
and kinetic energy of the impinging Na atom. We have found that the neutral Na 
is not likely to penetrate into the Ar matrix and sticks to the Ar surface for initial 
kinetic energy lower than ~ 0.2 meV while being reflected for larger impact kinetic 
energies. In case of the positively charged Na + , inclusion is observed, whatever 
the initial kinetic energy. The Ar matrix (or the finite Ar cluster) then undergoes 
strong perturbations and ejects one or more atoms to create a vacancy for the 
Na+ inclusion. As a first exploratory example, we have also studied the deposition 
of a neutral Nag cluster. As for the neutral atom, we see also a binding to the 
surface and no penetration into the Ar substrate. Somewhat surprisingly, there is 
only little perturbation of the Na cluster internal structure. Continued systematic 
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investigations on metal cluster deposition in Ar substrate are in progress. 
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Appendix A. The Na-Ar energy functional in detail 

The degrees of freedom of the model are: 

{<£ ra (r), n = l...A i} wavefunctions of cluster electrons 

{Rj, / = l...Mon} coordinates of cluster's Na+ ions 

{R a , a = l...A?Ar} coordinates of Ar atoms (cores Ar^ + ) 

{R' a , a = l...AA r } coordinates of the Ar valence clouds 
An Ar atom is described by two constituents with opposite charge, positive Ar 
core and negative Ar valence cloud, which allows a correct description of polarization 
dynamics. We associate a Gaussian charge charge distribution to both constituents 
having a width of the order of the 3p shell in Ar, in the spirit of. 32 The dynamical 
polarizability of the Na + ions is neglected and we treat them simply as charged 
point particles. 

The total energy of the system is composed as: 



-Etotal = -^Naclustcr + ^ Ar + -^coupl + -EVdW , 



The energy of the Na cluster i?Naciustor consists out of TDLDA (with SIC) for the 
electrons, MD for ions, and a coupling of both by soft, local pseudo-potentials, 
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see 27,28 for details. The Ar system and its coupling to the clusters is described by 
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The various contributions arc calibrated from independent sources, with a final 
fine tuning to the NaAr dimer modifying only the term W e iAr- The parameters are 
summarized in the table. The third column of the table indicates the source for the 
parameters. 

The (most important) polarization potentials are described by a valence electron 
cloud oscillating against the raregas core ion. Its parameters arc: q Al the effective 
charge of valence cloud, — q Ar m e \ the effective mass of valence cloud, k Al the 
restoring force for dipoles, and a Ar the width of the core and valence clouds. The 
q Al and k Ar are adjusted to reproduce the dynamical polarizability a_o(w) of the Ar 
atom at low frequencies, namely the static limit a£>(u; = 0) and the second derivative 
of a'Jj(d' = 0). The width a Ar is determined consistently such that the restoring 
force from the folded Coulomb force (for small displacements) reproduces the spring 
constant k Al . 
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a Ar =11.08 a 3 , 
w = 1.755 Ry 


WdAr 


A cl =0.47 , /3 e i=1.6941/a , r e i =2.2 a 


fit to NaAr 


T/ (core) 
K ArAr 


A Ar =1.367* 1(T 3 Ry , i? Ar =6.501 a 


fit to bulk Ar 


^ArNa 


/3 Na = 1-7624 a^ 1 , a Na = 1.815 a , A Na = 334.85 






C Na , 6 = 52.5 a 6 , , C Na , 8 = 1383 ag 


after 31 



The short range repulsion is provided by the various core potentials. For the Ar- 
Ar core interaction we employ a Lennard-Jones type potential with parameters re- 
producing binding properties of bulk Ar. The Na-Ar core potential is chosen accord- 
ing to, 31 within properly avoiding double counting of the the dipolc polarization- 
potential. 

The pseudo-potential W e i Ar for the electron- Ar core repulsion has been modeled 
according to the proposal of 32 with a final slight adjustment to the Na-Ar dimcr, 
data taken from from 36 and. 37 

The Van-der-Waals energy i?vaw is a correlation from the dipole excitation 
in the Ar atom coupled with a dipole excitation in the cluster. We exploit that 
WMic "C A_E Ar which simplifies the term to the variance of the dipole operator 
in the cluster, using again the regularized dipole operator f a corresponding to the 
smoothened Ar charge distributions. The full dipole variance is simplified in terms 
of the local variance. 
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